Anomalous elasticity in a disordered layered XY model 
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We investigate the efTects of layered quenched disorder on the behavior of planar magnets, superflu- 
ids, and superconductors by performing large-scale Monte-Carlo simulations of a three-dimensional 
randomly layered XY model. Our data provide numerical evidence for the recently predicted 
anomalously elastic (sliding) intermediate phase between the conventional high-temperature and 
low-temperature phases. In this intermediate phase, the spin-wave stiffness perpendicular to the 
layers vanishes in the thermodynamic limit while the stiffness parallel to the layers as well as the 
spontaneous magnetization are nonzero. In addition, the susceptibility displays unconventional 
finite-size scaling properties. We compare our Monte-Carlo results with the theoretical predic- 
tions, and we discuss possible experiments in ultracold atomic gases, layered superconductors and 
in nanostructures. 
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I. INTRODUCTION 

Extended defects can be found in a wide variety of 
condensed matter systems. For example, realistic mate- 
rials often contain one-dimensional and two-dimensional 
defects in the form of dislocation lines or grain bound- 
aries. Recent progress in nano-technology also allows re- 
searchers to custom-design artificial structures with anal- 
ogous properties. Extended defects in systems of ul- 
tracold atomic gases can be created by means of one- 
dimensional or two-dimensional disordered optical lat- 
tices. 

Extended defects are larger than the usual finite-size 
impurities and are thus harder to "average out." Conse- 
quently, they have a much greater influence on the ther- 
modynamic behavior of the system in question. This was 
first established on the example of the McCoy- Wu model, 
a two-dimensional disordered classical Ising model whose 
disorder is perfectly correlated in one of the two dimen- 
sions, i.e., it takes the form of parallel line defects. Mc- 
Coy and Wu [1-4] demonstrated in a series of papers that 
this model exhibits an exotic phase transition at which 
the magnetic susceptibility is infinite over an entire tem- 
perature range while the specific heat is smooth. Fisher 

0] later used a strong-disorder renormalization group 
to show that the critical point is of infinite-randomness 
type and is acc omp anied by strong (power-law) Griffiths 
singularities 10]. Ising models with plane defects, i.e., 
with perfect disorder correlations in two rather than one 
dimensions display even stronger disorder effects: instead 
of being sharp, the phase transition is smeared over a 
range of temperatures [ill, Ell • 

The effective dimensionality of the defects forms the 
basis of a classification 'I'f of phase transitions in 
systems with quenched disorder. Three classes can be 
distinguished: (i) If the defect dimensionality is below 
the lower critical dimension of the problem, the dis- 
ordered system has a conventional critical point with ex- 
ponentially weak Griffiths singularities, (ii) If the defect 
dimensionality is exactly equal to the lower critical di- 



mension, the critical point of the disordered system is 
of infinite-randomness type and accompanied by strong 
power-law Griffiths singularities, (iii) If the defects are 
above the lower critical dimension, individual regions can 
order independently, leading to a smearing of the global 
phase transition. 

The case of two-dimensional (planar) defects in sys- 
tems with XY symmetry is of particular conceptual and 
experimental importance. Theoretically, XY systems 
with perfect disorder correlations in two dimensions are 
right at the boundary between cases (ii) and (iii) in the 
above classification. True long-range order on individ- 
ual two-dimensional "slabs" is forbidden by the Mermin- 
Wagner theorem fTE']; however, these regions undergo 
a Kosterlitz-Thoulcss transition to a quasi long-range 
ordered phase [HI]. Experimentally, order parameters 
with XY symmetry occur not only in planar magnets 
but also in superconductors and superfluids. The fate 
of the XY phase transition with two-dimensional defects 
is thus of great interest for magnetic and superconduct- 
ing multilayers as well as ultracold atomic gases in one- 
dimensional disordered optical lattices. 

Recently, two simultaneous publications [l3, [13 inves- 
tigated this question theoretically. They predicted that 
the conventional high-temperature and low-temperature 
phases of a randomly layered XY system are separated by 
an anomalously elastic intermediate phase. In this exotic 
"sliding" phase which is part of the Griffiths region, the 
spin-wave (or superfluid) stiffness parallel to the layers 
is nonzero while the stiffness perpendicular to the layers 
vanishes. 

In this paper, we report the results of large-scale 
Monte-Carlo simulations of a three-dimensional ran- 
domly layered XY model which provide support for the 
phase transition scenario predicted in Refs. ^,[15. In 
particular, we give numerical evidence for existence of 
the anomalously elastic intermediate phase. Our paper 
in organized as follows. In Sec. [ill we define the randomly 
layered XY model. We briefly summarize the theoreti- 
cal predictions in Sec. IIIIl Section IIVI is devoted to the 
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FIG. 1. (Color online) Sketch of the layered XY model 
Layers of two different materials (represented by different in- 
teraction strengths j") are arranged in a random sequence. 

Monte-Carlo simulations and their results. We conclude 
in Sec. H 

II. RANDOMLY LAYERED XY MODEL 

In the following, we formulate the problem in the lan- 
guage of the XY ferromagnet. The results will apply to 
all phase transitions having 0(2) or U(l) order parame- 
ters, if expressed in terms of the appropriate variables. 

We consider a three-dimensional magnet consisting of 
a random sequence of layers made up of two different 
materials as shown in Fig. [T] Its Hamiltonian is a three- 
dimensional classical XY model defined on a lattice of 
perpendicular size L_l (in the z direction) and in-plane 
size L|| (in the x and y directions). It reads 

H — — j| (Sr'Sr+x + Sr-Sr-l-y)— J'^ Sr'Sr-l-z. (1) 
r r 

Here, Sr is a two-component unit vector on lattice site 
r, and x, y, and z are the unit vectors in the coordinate 
directions. The interactions within the layers, j| , and 
between the layers, J^, are both positive and indepen- 
dent random functions of the perpendicular coordinate z. 
For simplicity, we take all to be identical, = 

p^ . while the j| are drawn from a binary probability 
distribution 

P( Jll ) = (1 - c) <5( Jll -Ju)+ c5{j\\ - Ji) (2) 

with Ju > Ji- Here, c is the concentration of the "weak" 
layers while 1 — c is the concentration of the "strong" 
layers. 

Let us discuss the thermodynamics of the randomly 
layered XY model H]) qualitatively. If the temperature 
T is above the upper Griffiths temperature T„ (defined as 
the critical temperature of a hypothetical clean system 
having j| = J,j for all z), the model is in a conventional 



FIG. 2. (Color online) Schematic behavior of the order pa- 
rameter (spontaneous magnetization) m and the stiffnesses 
Psjl and ps,± vs. temperature T for any bounded disorder 
distribution. SD and SO denote the conventional strongly dis- 
ordered and ordered phases, respectively. The Griffiths region 
(bounded by Tu and Ti) consists of the "non-anomalous" (G) 
and the anomalous (AG) Griffiths phases. For an unbounded 
distribution, T„ — >• oo. 

paramagnetic phase (denoted "strongly disordered" in 
the phase diagram shown in Fig. [2|). Analogously, for 
temperatures below the lower Griffiths temperature T/ 
(defined as the critical temperature of a system having 
j| = Ji for all z), the model is in a conventional ferro- 
magnetic ("strongly ordered") phase. The most interest- 
ing temperature range is the Griffiths phase between T„ 
and Ti. Here, rare thick slabs (rare regions) of strong 
(j| = Jit) layers can show local ferromagnetism while 
the bulk is still nonmagnetic. Individual such slabs can- 
not be truly long-range ordered but they can develop 
quasi long-range order via individual Kosterlitz-Thouless 
phase transitions (l6j . The exotic behavior predicted 
in Refs. [13, Ell arises from the interplay between the 
quenched disorder and the Kosterlitz-Thouless physics 
of these strongly interacting slabs. 



III. OPTIMAL FLUCTUATION THEORY 

In this section, we summarize the results of the opti- 
mal fluctuation theory of Ref . [13] . A slab (rare region) 
consisting of Lnn consecutive strong layers exists with an 
exponentially small probability w{Lft.R) ^ (1 — c)^^^ = 
exp(— c-Lj^/e) with c = — ln(l — c). It undergoes a 
Kosterlitz-Thouless phase transition at a temperature 
TKriLufi) which can be estimated from finite-size scal- 
ing via Tu - Tkt ~ L'^r"- (Here, v « 0.6717 ^ is the 
correlation-length critical exponent of a clean 3D planar 
(XY) magnet.) Gonsequently, at a given temperature T, 
all rare regions of thickness Lflfl < Lc{T) ^ {T^ — T)~'^ 
are (locally) in the paramagnetic phase while those hav- 
ing thicknesses Lrr > Lc{T) are in the quasi- long-range 
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ordered phase. 

We first consider the behavior of a single rare region 
(slab) in the quasi-long-range ordered phase. Accord- 
ing to Kosterlitz-Thouless theory [l^, its spin correla- 
tion function C(r) decays as a non- universal power of 
the distance |r|, 

C(r) ^ |rr" (|r| ^ ^) . (3) 

The exponent rj takes the value 1/4 right at the 
Kosterlitz-Thouless transition and behaves as l/i_R_R in 
the limit of Lufj — >■ oo. Its thickness dependence can thus 
be modeled as 'ri{Ljiji) = jLc{T)/ Lnn. The power-law 
correlations (jS]) lead to a nonlinear magnetization-field 
curve 

m ^ i7"/(4-r,) ^ 

which implies that the magnetic susceptibility of a single 
slab in the quasi-long-range ordered phase is infinite. 

We now turn to thermodynamic observables of the full, 
randomly layered system in the Griffiths phase T„ > T > 
T/. The spin- wave stiffness ps is defined via the work 
required to twist the spins on two opposite boundaries of 
the sample by a relative angle d. In the limit of small 9 
and large system size, the free-energy density / depends 
on 9 as 

fi9)-fiO) = lps{9/Lf (5) 

which defines ps- As our randomly layered XY model is 
anisotropic, we need to distinguish the parallel spin-wave 
stiffness pi from the perpendicular spin- wave stiffness pjr . 
To determine pi, we apply twisted boundary conditions 
at X = and x — L|| and set L = L|| in ([S]) while the 
boundary conditions are applied at z = and z = L±_ 
(using L = L±) to find p^. 

When the twist is applied in x-direction, all layers in 
the sample have same boundary conditions. The total 
free energy cost due to the twist is thus simply given 
by a sum over all layers. As only slabs that are in the 
quasi-long-range ordered phase have a nonzero stiffness, 
the total parallel stiffness reads 

/>oo 

Ps.\\ dLB.B.w{LRR) Ps^b.r{Lr.b) ■ (6) 

Here, Ps,rr{Lb.r) is the (parallel) stiffness of a single 
slab of thickness £_r_r. It is related to the value of the 
exponent rj via rj = T / {2-kPs^rr). Because the proba- 
bility w{Lbr) decays exponentially with increasing Lbr, 
the integral is dominated by its lower bound. To leading 
exponential accuracy, the parallel stiffness is thus given 
by 

~ exp[-cic(r)] exp[-a (r„ - Ty] (7) 

with a a non-universal constant. This means, it is non- 
zero anywhere in the Griffiths phase and develops an ex- 
ponential tail towards (see Fig. [5]). 



When the twist 9 is applied between the bottom (z = 
0) and the top {z = L±) of the sample, the local twists 
between consecutive layers will vary from layer to layer. 
Minimizing the elastic free energy with respect to these 
local twists leads to 

Ps.i_ ^ (8) 

where are the effective couplings between the rare 
regions, and (. . .) is the average over the sample. Be- 
cause the spatial positions of the rare regions are ran- 
dom, the distribution of their nearest-neighbor distances 
is a Poisson distribution, P{R) = Rkt 'Sxp{—R/ Rkt), 
where Rkt w(,Lc)~^ ^ exp[cLc(T)] is the typical sep- 
aration. The effective interaction between neighboring 
rare regions decays exponentially, J^{R) ^ exp(— i?/^o), 
where is the bulk correlation length. We thus arrive 
at a power-law distribution 

^(Jeff) - {Jtnf'-' ■ (9) 

of the effective interactions. The Griffiths dynamical ex- 
ponent z = Rkt/Co takes the value c» at T„, and de- 
creases with decreasing temperature. Using this distri- 
bution, we find that the average {1/J^) diverges as long 
as z > 1. This implies that the perpendicular stiffness 
Ps.± vanishes in part of the Griffiths phase, viz., between 
Tu and a temperature at which z reaches the value 1. 
In this temperature region, the elastic free energy density 
displays an anomalous dependence on the system size, 
/(e) - /(O) ~ L];^"^ corresponding to ps,± ~ -^i"^- 
For T < Tg the average (l/J^g) converges, leading to a 
nonzero perpendicular stiffness. Glose to Tg, the perpen- 
dicular stiffness is expected to behave as ps.± ^ Tg — T 
(see Fig. 12). 

Other quantities can be found along the same lines . 
For example, the spontaneous magnetization is nonzero 
for all T < Tu and shows a double-exponential tail to- 
wards the nonmagnetic phase. Close to r„, it takes the 
asymptotic form 

ln(m) ^ - exp[a(T„ - T)-"] , [T ^ T„-) . (10) 

If a magnetic field H is applied at temperatures T <Tu, 
the magnetization-field curve takes the unusual form 

\nim) ^ -^\\n{H)\{Tu - , (if ^ 0) , (11) 

for magnetizations larger than the double-exponentially 
small spontaneous magnetization ([TU|. 

For the comparison between theory and Monte-Carlo 
simulations, finite-size effects are an important issue. As 
an example, we discuss the dependence of the magnetic 
susceptibility in the Griffiths phase on the in-plane sys- 
tem size When L|| is finite, the susceptibility of a 
single slab in the quasi-long-range ordered phase is no 
longer infinite. Its L||-dependence can be obtained from 
integrating the correlation function ^ to the upper cut- 
off L|| which results in X-R-r(-^Ii) ~ ■ Summing this 



over all rare regions yields the total susceptibility (per 
unit volume) as 

X-^Lfexp{-[5Lc(T)ln(i||)]i/2} . (12) 

IV. MONTE-CARLO SIMULATIONS 

In this section, we report the results of Monte-Carlo 
simulations of the randomly layered XY model ^ by 
means of the highly efficient Wolff cluster algorithm [2l| . 
To capture the physics of the rare regions, we have simu- 
lated large system sizes of up to L± = 800 and L|| — 100. 
In the binary distribution ^ for the in-plane interactions 

j|, we have chosen the values J„ = 1 and Ji = 0.25. 
All the simulations have been performed for an impurity 
concentration c = 0.8. With these parameter choices, the 
Griffiths phase ranges from Ti « 0.97 to Tu « 2.20. All 
data are averages over a large number (from 100 to 300) 
of disorder realizations. For each realization we have used 
100 Monte-Carlo (Wolff) sweeps for equilibration and an- 
other 100 sweeps for measurements. 

To test the anomalous elastic properties predicted in 
part of the Griffiths phase, we have computed the paral- 
lel and perpendicular spin-wave stiffnesses. Finding the 
stiffnesses by actually performing simulation runs with 
twisted boundary conditions is not very efficient. How- 
ever, the stiffnesses can be rewritten in terms of expecta- 
tion values calculated in a conventional run with periodic 
boundary conditions. In the case of the perpendicular 
stiffness, the resulting formula [l^l (see also Q) reads 

Ps,±-^ ^ Jr,r'(Sr-SrO(^-^')' (13) 
{r,r') 




where N — L±L'^^ is the total number of sites, and k is a 
unit vector perpendicular to the plane of the XY spins. 
For the calculation of ||, the term (z — z') needs to be 
replaced by {x — x'). 

Using this formula, we have calculated the parallel and 
perpendicular stiffnesses of a system of sizes L± = 800 
and L|| = 100. The results are shown in Fig. [3] The 
two stiffnesses clearly behave differently. The perpendic- 
ular stiffness yOs,_L vanishes at a temperature Tg « 1.15 
while the parallel stiffness p^ || remains nonzero to signif- 
icantly higher temperatures and develops a tail towards 
Tu- In agreement with the theoretical predictions, we 
thus find an intermediate anomalously elastic (sliding) 
phase in which ps.± — but || ^ 0. To further test 
the theory, we plot ln(ps ||) vs. (T„ — T)^" in the inset of 
Fig. [31 According to (O, the data sufficiently close to 
should fall onto a straight line. As the inset shows, our 
results follow the prediction over more than 1.5 orders 
of magnitude in || (down to the resolution limit set by 
the Monte-Carlo noise). 




FIG. 3. (Color online) Perpendicular and parallel spin-wave 
stiffnesses, ps,± and ps,\\, as functions of temperature T for 
a system with sizes L± — 800 and = 100. The data are 
averaged over 100 disorder configurations. The inset shows 
that Pa II follows the prediction (O for T ^ Tu- 
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[ln(L„)]"^ 

FIG. 4. (Color online) Susceptibility x as a function of in- 
plane system size Ly for several temperatures in the Griffiths 
region. The perpendicular size is L± = 800; the data are 
averages over 300 disorder configurations. The solid lines are 
fits to the theoretical prediction ([12]). 



In addition to the spin-wave stiffnesses, we have also 
analyzed the finite-size behavior of the magnetic suscep- 
tibility. Figure [Jj shows the susceptibility x as a function 
of L|| for several temperatures in the Griffiths region. We 
have used system sizes — 800 and Ly = 10 to 100 be- 
cause the condition L± ^ L\\ (such that L±^ is effectively 
infinite) needs to be fulfilled when studying the depen- 
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dence of x on L\\ . 

To compare the simulations to the theory, we plot the 
Monte-Carlo data as ln(x/L|) vs. [ln(L||)]^/^. In such 
a plot, the functional form ([T^ yields a straight line. 
Figure m shows that our data are in good agreement with 
the theoretical prediction over a wide temperature range. 
Some deviations appear for temperatures close to T„ and 
large L||. They are likely caused by the fact that the 
perpendicular size L±^ is not truly infinite in our simu- 
lations. Thus, a typical sample will not contain any of 
the very thick rare regions that dominate x ^-t temper- 
atures close to Tu- From fits of the data to (fT2|) . one 
can obtain estimates of the cutoff length scale Lc{T). Its 
temperature dependence does not agree very well with 
the prediction Lc ^ {Tu — T)~'^ , probably because the 
theory holds asymptotically close to T„ while the simula- 
tions for such temperatures suffer from the finite system 
size Lj^, as discussed above. 

V. CONCLUSIONS 

In summary, we have performed large-scale Monte- 
Carlo simulations of a classical three-dimensional XY 
model with layered randomness. Our results provide 
support for the recentl y p redicted unconventional phase 
transition scenario (TtI Il8l|. In particular, we have found 
evidence for the existence of an anomalously elastic (slid- 
ing) intermediate phase between T„ and Tg. In this 
phase, the stiffness parallel to the layered randomness 
is nonzero while the perpendicular stiffness vanishes. We 
have also confirmed the unusual finite-size scaling behav- 
ior of the magnetic susceptibility. 

It is interesting to compare the present results for a 
randomly layered XY model with corresponding results 
for Ising and Heisenberg spins. In a randomly layered 
Ising model, each rare region (slab) corresponds to a two- 
dimensional Ising model; it can thus undergo a transition 
to a long-range ordered state independently of the bulk 
system. The global phase transition is therefore smeared 
[111 [l3 l . (In the Ising case, the tail of the magnetization 
towards decays as a single exponential, i.e., m is signif- 
icantly larger than in the double exponential tail (llOp V 
In contrast, for Heisenberg symmetry, the rare regions 
correspond to two-dimensional Heisenberg models which 
cannot develop long-range order independently. As a re- 
sult, the global phase transition is sharp. However, the 
layered randomness leads to exotic infinite-randomness 
critical behavior [2^, [23] . In our case of the randomly 
layered XY model, the rare regions can undergo a phase 
transition independently from each other, but only to a 



quasi-long-range ordered state rather than to true long- 
range order. In terms of the classification of phase 
transitions in disordered systems, the randomly layered 
XY model is thus features a hybrid between a smeared 
and a sharp transition. 

Turning to experiment, our results are applicable to 
a variety of systems. Even though our theory is formu- 
lated in the language of the planar (XY) ferromagnet, it 
holds for all thermal phase transitions with 0(2) or U(l) 
order parameters, if it is expressed in terms of the appro- 
priate variables. For randomly layered superconductors 
and superfluids, the magnetization should be exchanged 
for the Cooper pair amplitude or the condensate wave 
function. Analogously, the spin-wave stiffness should be 
substituted by the superfluid density. It is worth noting 
that recent large-scale quantum Monte-Carlo simulations 
of the Bose-Hubbard Hamiltonian on randomly stacked 
layers [25j also confirmed the existence of an anomalously 
elastic (sliding) phase between the normal fiuid and the 
superfluid in this system. 

Experiments in ultracold atomic gases have al- 
ready demonstrated the Kosterlitz-Thouless transition 
in stacks of quasi two-dimensional layers created by a 
strong one-dimensional optical lattice [1^. Moreover, 
disordered optical lattices have been used to study An- 
derson localization of matter waves [2?, '28']. Our results 
apply to large irregular stacks of quasi two-dimensional 
layers created by a one-dimensional disordered optical 
lattice of a strength that still allows some weak coupling 
between the layers. 

Recently, experiments on several layered perovskite su- 
perconductors H^-Ull found unexpected anisotropics of 
the superconducting properties that imply an apparent 
decoupling of the superconducting layers. Our anoma- 
lously elastic phase may explain these observations pro- 
vided that there is sufficient c-axis disorder is the sam- 
ples. Our theory could also be tested by manufacturing 
layered nanostructures of different magnetic or supercon- 
ducting materials. Magnetic multilayers with systematic 
variations of the critical temperature from layer to layer 
have already been produced .32j. The system we have 
studied can be realized as a random version of such a 
structure using an easy-plane magnetic material. 
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